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Abstract 

Dendrites with developed sidebranches are numerically studied with 
a coupled map lattice model. The competitive dynamics among side- 
branches determines the shape of the envelope. The envelope has a 
parabolic shape near the tip of the dendrite and the envelope angle with 
respect to the main branch increases up to 45° finally. In an interme- 
diate region, the envelope grows roughly in a power law, however, the 
exponent increases gradually as a function of the distance from the tip. 
The competitive dynamics among many branches is also observed in a 
unidirectional growth from a linear seed, and it is compared with the 
competitive dynamics of sidebranches. 

1 Introduction 

Pattern formations in diffusion fields have been studied as a typical nonlinear 
and noncquilibrium problem. The dendrite is a typical pattern observed in 
crystal growth. The tip velocity and the shape near the tip region have been 
intensively studied. The growth law of the dendrite is expressed as vp 2 — const, 
where v is the tip velocity and p is the radius of curvature of the parabolic 
tip. 00 The parabolic structure of the dendritic tip is linearly stable and the 
formation of sidebranches is considered to be due to some noise effects. [5] Per- 
turbations are expected to grow as exp^a; 1 / 4 ), where x is the distance from the 
dendritic tip,|H] and they grow into sidebranches. Well-developed sidebranches 
can be easily observed in experiments of succinonitrile and NH4CI. There are 
some quantitative investigations about the well-developed sidebranches. Huang 
and Glicksman investigated the development of sidebranches in an experiment 
of succinonitrile. [7] They found that the sidebranching spacing, which is the 
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interval between neighboring active sidebranches, grows in a power law of ex- 
ponent 1.3 as a function of distance from the dendritic tip. They found the 
angle of the envelope with respect to the main branch increases as a function 
of the distance from the tip until 41°. They expected that the angle increases 
up to 45° finally where the sidebranches grow with the same velocity as the 
main branch. Li and Beckermann analyzed the shape of the envelope of the 
sidebranches for succinonitrile on the microgravity flight and found that the 
envelope grows with a power law of exponent 0.859 as a function of the dis- 
tance from the dendritic tip. [H] The exponent is less than 1, and it implies that 
the envelope angle decreases as a function of distance from the tip, which is 
different from the observation by Huang and Glicksman. Recently Honda and 
Honjo are investigating the shape of the sidebranches of kinetic dendrites in 
viscous fingering on the plate with many grooves of square symmetry. [9] Theo- 
retically well-developed sidebranches have not been studied in detail. We will 
study well-developed sidebranches using a coupled map lattice model. 

We proposed a coupled map lattice model as a simple simulation method 
for crystal growth and succeeded in generating various patterns such as DLA, 
DBM and dendrites. [TO] The dendrite with anisotropic surface tension obeys the 
growth law vp 2 — const also in our coupled map lattice, when the undercooling 
is sufficiently small. m| The model is a deterministic model similar to the phase 
field model. ^2] 0m coupled map lattice models involves the kinetic effect and 
we have not succeeded in removing the kinetic effect as in the phase-field model 
by Karma and Rappel.^l However, our coupled map model is much simpler and 
numerically preferable for a large scale simulation even at small undercooling 
such as A ~ 0.1. Previously, different types of coupled map lattice models were 
proposed by Kaneko and Kessler et al.^jE] Their coupled map lattice models 
may be interpreted as direct discretization of the phase-field models, but they 
have not reported dendritic growth. Our model is a simpler model, in that only 
the linear diffusion process is assumed in lattice points other than interface sites. 
Packard and Li-Godenfeld proposed cellular automaton models for the crystal 
growth. [TBI 116) In their cellular automaton models, the order parameter takes 
(liquid) or 1 (solid). In our coupled map lattice model, the order parameter takes 
a continuous value between and 1 at the interface sites. The interface grows 
therefore more smoothly than the cellular automaton models. Our coupled map 
lattice model is considered to be a preferable model compared to the previous 
coupled map lattice models and cellular automaton models. We will perform 
numerical simulation of dendrites using the coupled map lattice model to study 
the growth law of the envelope. 

2 Sidebranch structure of dendrites in a coupled 
map lattice 

We consider melt growth of crystal with four-fold rotational symmetry. Our 
model consists of two steps of time evolution on a square lattice. The four-fold 
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Figure 1: Dendritic patterns in a coupled map lattice for (a) C2 = 0.3 and (b) 
0.7. The dendritic patterns are plotted as the tip position locates at the origin. 



rotational symmetry of crystal is involved in performing the simulation on the 
square lattice. If the triangular lattice is used, a pattern like a snowflake with 
the six- fold rotational symmetry grows as shown in the previous paper. [TT] The 
first step is a diffusion process and the second step is a growth process at the 
interface. The diffusion process of the latent heat released at the interface is 
expressed as 

u'n+S-J) = u„(i,j)+D{u„(i + l,j)+u n (i-l,j) + 

u n (i,j + 1) + u n (i,j - 1) - 4u„(i,i)}, (1) 

where u, u' denotes a dimensionless temperature, n is the step number, de- 
notes the lattice point and D is a diffusion constant. This process is interpreted 
as the simplest discretization of the Laplacian. The second step is a growth 
process at the interface. The order parameter is expressed with the variable 
x(i,j), where x(i,j) is in the liquid phase and 1 in the solid phase. The order 
parameter x(i,j) changes only at the interface sites between the liquid sites and 
the solid sites. The growth rule is written as 

x n +l(i,j) = x n (i,j)+c 1 (l~u' n+1 (i,j)), 

u n +i{i,j) = u' n+1 (i,j)+c 2 (l-u' n+1 (i,j)), (2) 

where the melting temperature is assumed to be I. There are two parameters 
ci and C2 in our model, and the ratio c^jcx means the latent heat. This growth 
process is done only at the interface sites, and x n +i, u n +\ for the other sites keep 
the values x n ,u' n+1 at the previous step. When x(i,j) goes over a threshold 
x c = 1, the interface site changes into a solid site. The diffusion process (1) 
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Figure 2: (a) Envelope profiles for C2 = 0.3,0.5,0.7 and 0.9, which are plotted 
with the coordinates scaled with the diffusion length 2D /v. The numbers in the 
figure denote the values of ci (b) Double-logarithmic plots of the envelopes for 
C2 = 0.7 and 0.9. Two linear lines denote power laws with exponents 0.5 and 
1.4. 



and the growth process (2) at the interface are performed alternately. The 
initial temperature is assumed to be u(i,j) — uq = 0. Then, the normalized 
undercooling A is given by A = (1 — Uo)/ (02/01) — c\jci. We have set a seed of 
crystal at the origin as an initial condition. To save computation time, we have 
calculated only the upper half of the crystal assuming the mirror symmetry. 
The system size is 1200 x 1200 in most simulations. The parameter c\ is fixed 
to be 0.1 and the diffusion constant is D = 0.2. The surface tension effect is 
neglected in this paper. (A stable parabolic structure is obtained if the term for 
the surface tension is included in this model as shown in ref.ll.) Our dendrite 
is therefore a kinetic dendrite, since the growth direction is determined by the 
kinetic anisotropy.|17| Our simulation may be closely related to the experiment 
by Honda and Honjo, since they investigate kinetic dendrites in two dimensions. 

Figure 1 displays snapshot patterns of dendrites for (a) C2 = 0.3 and (b) 
C2 = 0.7. To include small noise effect, we have assumed that the threshold 
x c at the tip site (only one site) takes a random value between 0.99 and 1.01. 
Sidebranches are well developed in the region far from the tip. The sidebranches 
which locate in x > 1000 are not shown. The sidebranches are more developed 
for C2 = 0.3 than for C2 = 0.7, since the undercooling is large. We have calcu- 
lated the envelope line for each snapshot pattern and constructed an average 
curve of the envelope for c 2 = 0.3, 0.5, 0.7 and 0.9. In our numerical simulation, 
an envelope of sidebranches is defined as the line which connects the tip posi- 
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Figure 3: The derivative d\ny/d\n function of the scaled distance from 

the tip x for c 2 = 0.5, 0.7 and 0.9. 



tions of active sidebranches. A sidebranch is judged as active if it is larger than 
the neighboring active sidebranch locating in the left of itself (with smaller i). 
As a result of our definition, the envelope is a monotonously increasing func- 
tion. A relevant length scale of our model is the diffusion length I = 2D /v. The 
tip velocities are numerically obtained respectively as 0.0058,0.0081,0.0124 and 
0.0227 for c 2 = 0.9, 0.7, 0.5 and 0.3. The velocities do not obey the law v ~ A 4 
like the surface tension dendrite in two dimensions, but seem to obey v ~ A 
for the kinetic dendrites in our model. Figure 2(a) displays the averaged enve- 
lope, which is plotted with the scaled length unit x = i/l,y — j/L (We have 
performed numerical simulation with larger size 1500 x 1500 only for c 2 = 0.9, 
since the development of sidebranches is slow at the parameter.) The envelope 
for c 2 = 0.3 can be approximated as y — x — 14.5 for large x. That means 
that the angle of the envelope with respect to the x-direction approaches 45°. 
The envelopes overlap with each other rather well, especially, the envelope for 
c 2 = 0.9 overlaps with that for c 2 = 0.7 very well. There may exist a universal 
scaling function for the envelope function for sufficiently small undercooling. 
Figure 2(b) displays logarithmic plots of the envelope functions for c 2 — 0.7 and 
0.9. For small x, the envelope is approximated as x 1 / 2 , where the sidebranches 
are not well developed and the parabolic structure of the tip region is main- 
tained. For relatively larger x, the envelope, roughly speaking, obeys a power 
law with exponent about 1.4. However, the envelope in the logarithmic scale 
has a tendency to curve upwards as x is increased. The derivative d In y/d In x 
is calculated and displayed in Fig. 3 for c 2 = 0.5,0.7 and 0.9. The derivative 
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Figure 4: Sidebranch spacing d(x) as a function of the scaled distance from the 
tip x for c 2 = 0.9 (solid curve) and 0.7 (dashed curve). 

is about 0.5 for small x, which corresponds to the parabolic shape, however, 
it increases from 1.2 to 1.9 in the intermediate region as x is increased. This 
implies that the envelope does not obey a simple power law. We have also cal- 
culated the spacing between neighboring active sidebranches and obtained the 
average curve of the spacing between the active sidebranches as a function of 
the distance from the tip position. Figure 4 displays the average spacing d(x) 
for c 2 = 0.7 and 0.9 as a function of the scaled distance x. The two curves again 
overlap very well. The sidebranch spacing is almost a linear function of x. By 
the fitting to a power law d(x) ~ x 13 , the exponent (5 is evaluated to be about 
1.15. The sidebranch spacing increases with the distance from the tip position 
because of the competition among sidebranches. Smaller sidebranches stop to 
grow and the survived sidebranches become larger. Competitive processes be- 
tween survived large sidebranches go on further, until the spacing between the 
sidebranches becomes sufficiently larger than the diffusion length. 

3 Competitive dynamics among dendritic branches 
in the unidirectional growth from a linear seed 

We can see the competitive growth of dendritic sidebranches more explicitly 
using another type simulation, which is the crystal growth starting from a linear 
seed. We have performed a numerical simulation in a square box of L x x L y — 
600 x 600. We have imposed periodic boundary conditions in the i direction such 
as u(L xl j) — u(l,j),u(0,j) = u(L x — at i — L x and 0, and fixed boundary 
conditions such as u(i,L y ) — at j — L y and no-flux boundary conditions 
u(i,0) = u(i, 1) at j — 0. As an initial condition for a linear seed, we assumed 
that x{i, 1) = l,u(i, 1) = 1, x{i,2) = 0, and u{i,2) took a random number 
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Figure 5: (a) Snapshot pattern of dendritic branches for the unidirectional 
growth from a linear seed for C2 = 0.7. (b) Contour curves of u = 
0.1,0.5,0.9,0.95,0.99 and 0.999 for the temperature field u at the same time 
as (a). 



between and 0.1. The initial value of x(i,j) and u(i,j) for j > 3 are all 0. 
Figure 5(a) is a snapshot pattern of x(i,j) at c 2 = 0.7. Many small dendritic 
branches are created initially, but the number of growing branches decreases 
in time by the competitive interaction between the neighboring branches. The 
velocities of delayed branches decrease and finally become almost zero. Large 
branches seen in Fig. 5(a) are survived branches for the competitive process. 
Figure 5(b) displays contour lines of u = 0.1, 0.5, 0.9, 0.95, 0.99 and 0.999. The 
large neighboring branches interact competitively through the diffusion field u. 
The temperature u is almost 1 in the regions between the delayed branches. 

Figure 6(a) displays the time evolution of the tip postition j p of the largest 
branch as a function of the step number n. The ensemble average with respect 
to 20 different initial conditions is taken. The form of this time evolution is 
similar to the envelope of sidebranches shown in Fig. 2. We have compared the 
two curves by rescaling the coordinates. The step number n can be converted 
into a length by i = vn where v is the tip velocity of the main stem and the 
scaled coordinates x = i/(2D/v) and y = j/(2D/v) are used in Fig. 6(b). 
Figure 6(b) displays the relation of y — 0.67 j p / (2Dv) and x — nvj (2D/v) (solid 
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Figure 6: (a) Time evolution of the tip position j p of the largest branch for 
c 2 = 0.7. (b) Comparison of the time evolution of the tip position of the largest 
branch for the unidirectional growth from a linear seed (solid curve) with the 
envelope shown in Fig. 2 for c 2 = 0.7 (dashed curve). They overlap very well.(c) 
Comparison of the time evolution of the spacing between the neighboring large 
branches for the unidirectional growth (solid curve) with the spacing between 
the active sidebranches shown in Fig. 4 for c 2 = 0.7 (dashed curve). 

curve) and it is compared with the envelope of sidebranches at c 2 = 0.7 (dashed 
curve). The two plots overlap very well, although the numerical factor 0.67 
is suitably chosen. We have also calculated the time evolution of the average 
interval d(n) between active branches for the crystal growth starting from a 
linear seed. Here, the active branches at step n are chosen as the branches whose 
tip positions are larger than 0.7 j p . The average interval between active branches 
increases almost in proportion to n, which corresponds to the linear growth of 
the sidebranch spacing shown in Fig. 4. Figure 6(c) compares the time evolution 
of the rescaled quantity d(x) = 0.8d(n)/(2D/v) as a function of x = nvj (2D/v) 
(solid curve) with the average spacing between active sidebranches shown in 
Fig. 4 for c 2 = 0.7 (dashed curve). (The numerical factor 0.8 is again suitably 
chosen.) The good agreement of two plots in Fig. 6(b) and 6(c) implies that the 
competitive dynamics among dendritic branches starting from a linear seed and 
that for sidebranches in the region far from the dendritic tip are essentially the 
same. 

The tip velocity of the largest branch is small at small n and it increases with 
time as is seen in Fig. 6(a). This is because the spacings between the survived 
branches increase with time and the tip velocities become larger for the larger 
spacing. To understand the relation between the tip velocity and the spacing 
between the branches, we have performed numerical simulations in a rectangular 
box of size L x x L y with periodic boundary conditions in the direction of x. As an 
initial condition, we have put a small branch of size 3 x 10 at i = L x /2 and j = 1. 
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Figure 7: (a) Time evolution of the tip position of one branch in a rectangular 
box with size 50 x 700 for = 0.7. The linear line denotes the power law with 
exponent 0.7. (b) Numerically evaluated exponent a of the power-law growth 
as a function of the width L x of the rectangular box. 

Since we have assumed periodic boundary conditions, the growth condition is 
equivalent to the case involving many branches, which have completely the same 
size and are periodically located with the same interval L x , in an infinitely wide 
system. The system size L y is fixed to be 700, and L x is changed from 10 to 120, 
and the parameter c-i is 0.7. Figure 7(a) displays the time evolution of the tip 
position of the branch in a logarithmic scale. The time evolution seems to obey 
a power law j p = n a with exponent 0.7 in fairly wide region of t. In a narrow 
channel, the steady growth of dendrite is not possible when the undercooling 
is small. Brener et al. investigated the parameter region for the existence 
of steadily growing dendrites. JS| Figure 7(b) displays numerically evaluated 
exponents a as a function of L x . The flat interface grows with a power law of 
exponent 1/2, which is shown theoretically and was numerically shown also in 
our coupled map lattice in the previous paper. [TT] which corresponds to the limit 
of L x — 0. Inversely, the exponent becomes 1, when the undercooling and the 
spacing L x are larger than a certain threshold, and the steady growth becomes 
possible. Figure 7(b) shows that the exponent increases with the spacing L x for 
our kinetic dendrites. 

In the crystal growth from a linear seed, the interface is almost flat and the 
growth exponent is nearly 1 /2 in an initial stage. Then, many small branches are 
created. At this stage, the branch spacing is small, the growth is slow, and the 
growth exponent is still close to 1/2. As the competitive process proceeds, the 
average spacing between active branches increases as a result of competition, 
and the tip velocities of the survived branches increase further. In this time 
scale, the tip position seems to grow like a power law with exponent larger than 
1. The exponents of the power law may be evaluated roughly as follows. If the 
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tip velocity vp grows in a power law with exponents a as shown in Fig. 7(b), the 
exponents depend on the average spacing I, and the average spacing I increases 
roughly as f ~ t, the tip position will grow as jp <~ J vpdt <~ J d/dt{t a ^}dt ~ 
J dt{at a ~ 1 + (da/dt)t a \nt}. If a grows as a ~ 0.5 + ci for small i with a certain 
constant c, jp ~ t and the exponent for the tip position is evaluated as 1.5 
for small t, although this is a very rough approximation. If the spacing becomes 
comparable or larger than the diffusion length, the tip velocities of survived 
branches approach the constant tip velocity for a single branch in an infinitely 
wide system. 

4 Summary 

We have performed a large scale numerical simulation of kinetic dendrites with 
a coupled map lattice model. Well-developed sidebranches have appeared in 
the coupled map lattice. Many sidebranches are naturally created and the 
competitive dynamics among sidebranches determines the shape of the envelope. 
The envelope has a parabolic shape near the tip of the dendrite, and the angle 
of the envelope with respect to the main stem grows up to 45° finally. In 
an intermediate region, the envelope grows, roughly speaking, in a power law, 
however, the exponent seems to increase gradually as a function of the distance 
from the tip. 

Similar competitive dynamics among many dendritic branches is observed 
in the unidirectional growth. We could compare quantitatively the temporal 
evolution of the unidirectional growth with the envelope of sidebranches as a 
function of the distance from the tip position. The temporal evolution of a single 
branch in a narrow rectangular box seems to obey a power law in a wide range 
of t. The exponent increases with the width of the rectangular box. The width 
can be interpreted to correspond to the spacing betweeen the sidebranches for 
the dendritic growth. As a result of the competition among sidebranches, the 
spacing between the active sidebranches increases, and then the tip velocities of 
the sidebranches are accelerated further. When the spacing becomes sufficiently 
large, the tip velocities of the active branches become nearly equal to the tip 
velocity of a single branch in an infinitely wide system. 

We have calculated numerically the spacing d(x) and the envelope y(x) of 
sidebranches. The competitive dynamics among the active sidebranches via 
the diffusion field determines d(x) and y(x). Expecting that the large scale 
competitive dynamics may not depend on the surface tension effect and we 
have neglected it in this paper. The surface tension is important to determine 
the tip radius and the tip velocity as in the experiments by Glicksman et al. 
It is not difficult to involve the surface tension effect in our model. However, 
the surface tension effect makes the parabolic shape of the dendrites stable, 
and strong noises are necessary to construct well-developed sidebranches. It 
is a future problem to study the sidebranch structure for the surface tension 
dendrites. 
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